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We critically review several recent approaches to solving the two cosmo- 
logical constant problems. The "old" problem is the discrepancy between the 
observed value of p\ and the large values suggested by particle physics models. 
The second problem is the "time coincidence" between the epoch of galaxy 
formation to and the epoch of A-domination t&. It is conceivable that the 
"old" problem can be resolved by fundamental physics alone, but we argue 
that in order to explain the " time coincidence" we must account for anthropic 
selection effects. Our main focus here is on the discrete-A models in which 
A can change through nucleation of branes. We consider the cosmology of 
this type of models in the context of inflation and discuss the observational 
constraints on the model parameters. The issue of multiple brane nucleation 
raised by Feng et. al. is discussed in some detail. We also review continuous- 
A models in which the role of the cosmological constant is played by a slowly 
varying potential of a scalar field. We find that both continuous and discrete 
models can in principle solve both cosmological constant problems, although 
the required values of the parameters do not appear very natural. M-theory- 
motivated brane models, in which the brane tension is determined by the 
brane coupling to the four-form field, do not seem to be viable, except per- 
haps in a very tight corner of the parameter space. Finally, we point out that 
the time coincidence can also be explained in models where A is fixed, but 
the primordial density contrast Q = 5p/p is treated as a random variable. 



I. INTRODUCTION 

The cosmological constant A presents us with at least two intriguing problems. Particle 
physics models suggest that the natural value for this constant is set by the Planck scale, 
Mp ~ 10 18 GeV [we use the reduced Planck mass Mp = (87rG) -1 / 2 ]. The corresponding 
vacuum energy density is pa ~ Mp, which is some 120 orders of magnitude greater than the 
observational bounds. In supersymmetric theories, one can expect a lower value, 

Pa ~ Vsusy, (!) 

where t/susy is the supersymmetry breaking scale. However, with rjsusY ~ 1 TeV, this is still 
60 orders of magnitude too high. Until recently, this discrepancy between the expected and 
observed values was the only cosmological constant problem. Its solution, many believed, was 
that something so small could only be zero, due to some unknown symmetry or dynamical 
cancellation. 
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Thus, it came as a surprise when recent observations [IJ provided evidence that the uni- 
verse is accelerating, rather than decelerating, suggesting a non-zero cosmological constant^. 
The observationally suggested values of A correspond to p\ ~ pmo, where pmo is the present 
density of matter. This brings yet another puzzle. It is difficult to understand why we hap- 
pen to live at the epoch when pu ~ Pa- Another statement of the problem is why the time 
when A starts dominating the universe nearly coincides with the epoch of galaxy formation, 

t A ~ t G . (2) 

This is the so-called cosmic coincidence problem. 

A number of proposed solutions to these problems have recently appeared in the literature 
1^-0. Some of them rely on some form of the anthropic principle, while others do not. To 
our knowlege, the only approach that can explain both puzzles is the one that attributes 
them to anthropic selection effects. In this approach, the cosmological constant is assumed 
to be a random variable that can take different values in defferent parts of the universe. 

The purpose of this paper is to give a critical analysis of the proposed approaches, both 
anthropic and otherwise. Our main focus will be on the models with a discrete spectrum of A 
which have recently attracted much attention. We shall consider these models in the frame- 
work of inflationary cosmology and discuss the calculation of the probability distribution for 
Pa, as well as the observational constraints on the model parameters. 

The paper is organized as follows. In Section II we review the motivation for considering 
A as a random variable. In Section III we discuss models where A is a discrete variable, in 
particular the models where there is a four-form contribution to the cosmological constant, 
which may relax to a small value through nucleation of branes. In Section IV we analyze 
the cosmology of such models. In Section V we consider the possibility of coincident brane 
nucleation. In Section VI we discuss models where the cosmological constant is a continuous 
variable. In Section VII we consider the possibility of a slowly varying four-form field in 
theories with extra dimensions. In Section VIII we review some non-anthropic approaches 
to the problem. In Section IX we consider models where the time coincidence is explained 
by assuming that the primordial density contrast Q = 5p/p (and not necessarily A) is a 
random variable. Our conclusions are summarized in Section X. 



II. A AS A RANDOM VARIABLE 



Not all values of A are consistent with the existence of conscious observers. This observa- 
tion was made by Barrow and Tipler pi (see also [pH), but the first quantitative analysis is 
due to Weinberg |12|. In a spatially flat universe with a cosmological constant, gravitational 



l The surprise, however, was not total. In Ref. ^ (well before the supernova data JTJ would give 
the first observational evidence in this direction) it was noted that anthropic selection effects would 
place the cosmological constant in the range pa/ Pmo ~ 10, and that "the actual value is likely to 
be comparable to this upper bound." For a flat universe this implies 0^ ~ 0.9, not far from the 
observed value and certainly compatible with it, within the accuracy of the prediction. Similar 
predictions where made in S at about the same time. 
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clustering effectively stops at t ~ t A . At later times, the vacuum energy dominates and the 
universe enters a de Sitter stage of exponential expansion. An anthropic bound on pa can be 
obtained by requiring that it does not dominate before the redshift z max when the earliest 
galaxies are formed. Weinberg took z max ~ 4 and obtained 

p A £ 100pA/o- (3) 

This is a dramatic improvement over Eq. ([!]), but it still falls short of the observational bound 
by a factor of about 30. 

The anthropic bound @ specifies the value of pa which makes galaxy formation barely 
possible. However, as it was pointed out in @||, the observers are where the galaxies are, 
and thus most of the observers will detect not these marginal values, but rather the values 
that maximize the number of galaxies. More precisely, the probability distribution for pa 
can be written as 

dV( PA ) = V*{p A )v{p A )dp A . (4) 

Here, V*(p A )dp A is the a priori distribution, which is proportional to the volume of those 
parts of the universe where pa takes values in the interval dp A , and z^(pa) is the average 
number of galaxies that form per unit volume with a given value of pa- According to the 
"principle of mediocrity", which assumes that we are typical observers, Eq. (|4]) gives the 
probability distribution for us to observe a given value of pa- The calculation of v(ph) is 
a standard astrophysical problem; it can be done, for example, using the Press-Schechter 
formalism ||13|| . The a priori distribution P*(pa) should be determined from the theory of 
initial conditions, e.g., from an inflationary model. 

Martel, Shapiro and Weinberg |14| (see also fll5"l ) presented a detailed calculation of 



dV{p A ) assuming a flat a priori distribution, 

V*(p A ) = const (5) 

in the range of interest (|f). They found that the peak of the resulting probability distribution 
is close to the observationally suggested values of pa- The cosmic time coincidence is easy to 
understand in this approach |16|,17| if one notes that regions of the universe where t A ^ to 



do not form any galaxies at all, whereas regious where t A ^> tc are suppressed by "phase 
space" , since they correspond to a very tiny range of A. It was shown in Ref. that the 
probability distribution for tc/t A is peaked at ta/t A « 1.5, and thus most observers will find 
themselves in galaxies formed at to ~ t A . 

This anthropic solution to the cosmological constant problems is incomplete without a 
particle physics model that would allow A to take different values and without a theory 
of initial conditions, such as an inflationary cosmological model, that would allow one to 
calculate the a priori distribution V*(p A ). 

One possibility is to consider models in which the role of the vacuum energy is played 
by a slowly varying potential V(4>) of some scalar field 0, which is very weakly coupled to 
ordinary matter. The values of <fi are randomized by quantum fluctuations during inflation, 
and analysis shows that the resulting a priori distribution is indeed flat for a wide class 
of potentials PJl8|. The main challenge one has to face in this approach is to justify the 
exceedingly flat potential V(<j)) required by the model. We shall comment on this issue 
in Section VI. Before that, we shall consider an alternative possibility which has recently 
attracted much attention. This is provided by models with a discrete spectrum of pa- 
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III. MODELS WITH A DISCRETE SPECTRUM OF A 



The first model of this type was suggested in an early paper by Abbott |0| as an attempt 
to solve the old cosmological constant problem. He considered a self-interacting scalar field 
4> with a "washboard" potential V((p) of the form illustrated in Fig.l. The potential has 
local minima at <p n = nt] with n = 0, ±1, ±2, ... , separated from one another by barriers. 
The vacuum at = <p n has energy density 

P\n = ne + const (6) 

and can decay through bubble nucleation to the vacuum at <p n -i- 



V(0) 




FIG. 1. The washboard potential. 
The nucleation rate r„j per unit spacetime volume is given by 

T nl = A n e- B ", (7) 



where B n is the action of the Coleman - de Luccia instanton pC§ and the meaning of the 
subscript j will become clear shortly. The bubble radius at nucleation R n is bounded by 



< R n < H' 1 , where 

is the square of the expansion rate of de Sitter space filled with the vacuum <p n . The horizon 
radius and the curvature radius of that space are both equal to H~ x . 



4 



An analytic expression for B n can be given in the thin wall approximation, when 5 -C R n 
2uj . The general expression is somewhat cumbersome and we shall only consider the limiting 



cases of interest. For R n <C H n 1 , B n is given by the flat space expression ET 



approximately independent of n. In this regime R n ~ 3<r/e, so we should have aH n /e <C 1. 
In the opposite limit, aH n /e ^> 1, we have (H~ l — R n ) <C H' 1 and 

B (waii) w 2 7i 2 aH- 3 . (10) 

The vacuum energy difference e is unimportant in this case, and the action coincides with 
that for domain wall nucleation p2 |. The prefactor in Eq.(|7|) can be estimated as (see e.g. 



A n ~ d 2 ^. (11) 

Eqs.(P)-([Tl|) apply under the condition that the gravitational effect of the wall is negligible, 

a < M|# n . (12) 

Upward quantum jumps from n _ x to n are also possible pi| . The corresponding 
nucleation rate is 



r( n -i)T = exp 
For e <^ p n this can be approximated as 



24tt 2 M 4 ' J 



P \Pn Pn-1 



57T 2 e 



rv (13) 



r (n _ 1)T = exp y — —j^ j r ni , (14) 

where we have used Eq.(|) for H n . 

In order for the anthropic explanation to work, one needs 

e < Pmo ~ (10- 3 eV) 4 , (15) 

and in order to have successful baryogenesis, the energy density during inflation should 
exceed (1 TeV) 4 , which corresponds to 

H > 1(T 3 eV. (16) 

Combining this with Eq.(|14]), we see that the probabilities of upward and downward jumps 
in pa during inflation are nearly equal, except perhaps in the borderline case when 

H ~ e 1/4 ~ 1CT 3 eV. (17) 



An alternative discrete model, first discussed by Brown and Teitelboim pqj , assumes 
that the cosmological constant is due to a four-form field, 
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which can change its value through the nucleation of branes. The total vacuum energy 
density is given by 

PA = Pbare + F* /2, (19) 

where pbare < is the 'bare' cosmological constant at F — 0. The change of the field strength 
across the brane is 

AF = ±q, (20) 

where q = const is fixed by the model. The four-form model has recently attracted much 
attention because four-form fields with appropriate couplings to branes naturally arise in 
the context of M-theory. In this case the brane tension is |5|||] 



a = qMp/V2, (21) 

and the effective thickness of the branes is 5 ~ Mp 1 , so that the thin wall approximation is 
justified. 

At present we should have |F| ~ (— 2p bare ) l l 2 ) so that the bare cosmological constant is 
almost neutralized. Then, in the range of interest, the spectrum of pa is nearly equidistant, 
with a separation 

Ap A = e « (-2p bare ) 1/2 q, (22) 
and the model is very similar to the Abbott's "washboard" model. We expect 

\ Pbare \ > {\TeV)\ (23) 
and it follows from Eq.([l5|) that q £ 1CT 90 Mj, and 

o < (1(T 3 eV)\ (24) 



where we have used the relation (pl[) between a and q. Such small values of q and a 
may appear problematic, but in a recent paper [[J Feng, March-Russell, Sethi and Wilczek 
(FMSW) have argued that they can naturally arise due to non-perturbative effects in M- 
theory. With a and H satisfying the bounds ( |2~4"D and (0), the condition of negligible brane 
gravity (|T2|) is also satisfied, and thus Eqs.(PD, (^0|) can be used. 

With the aid of Eqs. (|2l1) -(p4|) it can be easily verified that the flat space bounce action 
(|) is bounded by § 

B [flat) < 1Q2> ( 25 ) 

This inequality is saturated for p bare ~ (1 TeV) 4 and 

e V4 ^ 0-V3 _ 10 -3 eV> (26) 
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If a and e significantly differ from these borderline values, then 5^,1 and brane nucleation is 
unsuppressed. A similar bound is obtained for the wall nucleation action (|10|) using Eqs.(P4|) 
and QTD, 

B (wM) < 2Q ( 27 j 

Here, the inequality is saturated for 

H ~ a 1/3 ~ 10~ 3 eV. (28) 

We note that Eqs.(pT|),(|2"2"|) apply only to models based on M-theory, and therefore the 
constraints (H,©, and (|27[) are also limited to this class of models. 



A different version of the four-form model has been developed by Bousso and Polchinski 
(BP) ||. They assume that several four-forms Fj are present, so that Eq . (|19D is replaced by 

PK = Pbare + \Y, F l ( 29 ) 
/ i 

The corresponding "charges" g« are not assumed to be very small, but BP have shown that 
with multiple four-forms the spectrum of the allowed values of p\ can be sufficiently dense 
to satisfy the condition ([T^) in the range of interest. However, the situation here is quite 
different from that in the FMSW model. As pointed out by the authors themselves, and 
further emphasized by Banks, Dine and Motl JHJ, the vacua with nearby values of pa have 
very different values of Fi and are expected to have very different physical properties. There 
is no reason to expect the a priori probabilities for these vacua to be similar. Moreover, 
the low energy physics in different vacua is likely to be different, so the process of galaxy 
formation and the types of life that can evolve will also differ. It appears therefore that the 
anthropic approach to solving the cosmological constant problems cannot be applied to this 
case p6] . 



IV. COSMOLOGY OF THE FOUR-FORM MODELS 

A. A priori distribution 

We shall now discuss the four-form models with brane nucleation in the context of in- 
flationary cosmology. The energy density of the universe during inflation can be expressed 
as 

Pn(x)=U(x)+PAn, (30) 

where U(x) is the potential of the inflaton field \i pAn is the cosmological constant contri- 
bution (|T9|) , and index n labels the vacuum energies corresponding to different values of the 
four-form field F. [The inflaton potential is generally F-dependent and has different forms 
U n (x) in defferent vacua |5,26|. Here we shall disregard this difference, assuming that the 



variation of U(x) is negligible in the narrow anthropic range of pa that will be of interest to 
us.] The minimum of U(x) is assumed to be at U m i n = 0. The spacetime during inflation is 
locally approximately de Sitter, 
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ds * = dt 2 _ e 2H n t dx 2^ 



(31) 



with H n (x) given by Eq.(§). 

A remarkable feature of inflation, which will play an important role in our discussion here, 
is that generically inflation never ends completely in the entire universe. The evolution of 
the inflaton field x is influenced by quantum fluctuations, and as a result thermalization 
does not occur simultaneously in different parts of the universe. In most of the models, one 
finds that at any time there are parts of the universe that are still inflating and that the 
total volume of inflating regions is growing with time P7|]28[| . This picture is often referred 
to as stochastic, or eternal, inflation. 




FIG. 2. Inflaton potential 

The full range of the field x can be divided into the "diffusion" , slow-roll, and thermal- 
ization parts, as illustrated in Fig. 2. In the diffusion range, x ~ Xq-> the inflaton dynamics 
is dominated by quantum fluctuations. It is this regime that is responsible for the eternal 
nature of inflation. In the slow-roll regime, Xq ~ X ~ X*i the inflaton rolls down its poten- 
tial. As it reaches the thermalization point x*i ft starts oscillating about the minimum of 
the potential, and its energy gets thermalized. The hypersurfaces x = X* are therefore the 
boundaries between inflating and thermalized regions of spacetime. These surfaces play the 
role of the big bang for the corresponding thermalized regions. There is typically an infinite 
number of such surfaces, each of them having an infinite volume. (For a discussion of the 



spacetime structure of inflationary universe see, e.g., p9|.) 

As the inflaton x fluctuates, rolls down its potential, and eventually thermalizes its 
energy, spherical branes nucleate at the rates ©,(OD changing the local values of the four- 
form field F. All possible values of p\ n will be taken on each infinite thermalization surface 
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£*, and the a priori probability V* n can be defined as the fraction of the volume of £* 
occupied by regions with vacuum energy density p\ n . [The volume fraction on an infinite 
hypersurface can be defined by calculating this fraction in a sphere of geodesic radius R and 
taking the limit R — > oo.] 

Brane nucleation can both decrease and increase the value of p\ n \ the corresponding 
nucleation rates are related by Eq.(|l3|). For A-lowering events, the bubble radius is initially 
smaller than the horizon H~ x and then grows in the comoving coordinates, while for A-raising 
events the radius is initially larger than the horizon and then decreases in the comoving 
coordinates. In both cases, with an appropriate definition of the nucleation time, the radius 
of the bubble nucleated at t — asymptotically approaches H~ 1 e Hnt |50fl . This means that 
the region affected by each nucleation event is a sphere of initial radius H~ x . For a comoving 
observer in vacuum n, the probabilities per unit time to witness a A-raising or lowering event 
are 

An 

^ = V n} —H~\ (32) 

47T 

K nl = T nl —H-\ (33) 
It follows from (|T3|) that these probabilities are related by 

«(n-l)f = «ni(/n-l//n)> ( 34 ) 

where 

f n = H~ 3 exp( — -H. (35) 

Consider an ensemble of comoving observers and let p n {t) be the fraction of observers 
in the n-th vacuum, where t is the proper time along the observers' world lines. The time 
evolution of p n is described by the equations 

dp n /dt = -(« nT + K ni )p n + K( n -i)tPn-l + K (n+l)iPn+l- (36) 

Let us assume for a moment that the inflaton potential remains unchanged, 

U(x) = const, (37) 

so that K n i and K n [ do not change with time. Then the solutions of (^) approach the 
stationary distribution 

p n oc oc ^exp V (38) 

We shall be mostly interested in the probability distribution in the anthropic range ([3|), 
where p\ n can be approximated by (§) with e from Eq.(^), and Eq.(^8|) takes the form 

( 87r2e \ 

p n oc exp I ~-^ n I • ( 39 ) 
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If inflation is well above the electroweak scale, H ^> 10 3 eV, then the distribution fl39|) is 
nearly flat in the anthropic range, 

p n ~ const. (40) 

The assumption (|37|) may or may not be a good approximation, depending on the shape 
of the potential U(x)- A simple example of a model where this approximation is adequate 
is a "new inflation" type model with a very flat potential in the diffusion range near the 
maximum of U(x) and a relatively steep decline to the minimum in the slow roll range. The 
distribution fl38| ) is established during the very long diffusion period, and then it does not 
change much during the slow roll period if the duration of the slow roll is shorter than the 
characteristic bubble nucleation time. Here we shall assume that the approximation (|37|) is 
justified. 

Can the distribution (^) be identified with the a priori probability distribution P* n ? 
The answer is "Yes, but only in a restricted class of models". An ensemble of comoving 
observers gives a comoving- volume distribution for p^n, which does not account for the fact 
that regions with different values of p\ n expand at different rates. The condition for this 
effect to be negligible is that brane nucleations should reshuffle the values of p\ n between 
different regions on a timescale tb which is much shorter than the time r# it takes for the 
differential expansion rate to significantly modify the distribution, 

tb < r H . (41) 

As we noted in Section II, the probabilities of upward and downward jumps in p\ should 
be nearly equal, except perhaps in the borderline case (|iTp. This means that the evolution 
of pa can be pictured as a random walk with steps taken on a timescale r ~ k^ 1 ~ H^T^ 1 . 
The anthropic range (S) comprises iV ~ 10 2 p M0 /e steps, and thus 



t b ~ N'H'T- 1 ~ 10" ( -4L? B 



) ^wf- < 42 > 

where we have used Eqs.(^) and ([0]). [In this discussion we have dropped the subscripts ] 
and I, since the upward and downward nucleation rates are nearly equal, and the subscript 
n since H n is nearly constant in the anthropic range.] 

The variation of the expansion rate in the anthropic range of p\ is 

Ne 

SH ~ — ^— , (43) 

and the time th can be estimated as 

t h ~ 1/5H ~ 10- 2 M P H/p M0 . (44) 
The condition r B r H can now be expressed as 

e < 10 — ' (45) 

Pmo 

Parameter values satisfying this condition can be readily found. 
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What happens in the opposite limit, when tb 3> t#? In this case the differential ex- 
pansion is important and the probabilities for faster expanding regions with higher values 
of pa are strongly enhanced. The predicted values of p\ should therefore be significantly 



higher than those obtained with a flat a priori distribution. Martel et.al. |14j have found 
that in the latter case the probability distribution is peaked at p\ = 0, the width of the peak 
being somewhat broader than the observationally suggested value. Models with r# t# 
will have the peak displaced towards higher values of pa and are therefore unlikely to give 
a good agreement with observations. A quantitative analysis of probability distributions 
in such models can be given by a relatively straightforward generalization of the formalism 
developed in Ref. pjf . 

We note finally that in models with borderline values of parameters ( |17|) the a priori 
distribution ( p9|) can significantly deviate from flatness, with smaller values of pa being 
favored. This would displace the peak of the resulting distribution to negative values of p\ 
and if anything would make the observational situation even worse. 



B. Observational constraints 

Models of the type we are discussing here suggest that we live in a bubble surrounded by 
an expanding brane. The values of pa inside and outside the brane are different. Let us first 
assume that the visible universe is contained within a single bubble. This means that the 
brane surrounding our bubble nucleated before the presently observable universe crossed the 
horizon during inflation. For this situation to be typical, the brane nucleation rates should 
be rather low, both during inflation and at present. This requires that the corresponding 
bounce actions should be large, B ^> 1. In M-theory motivated models this is possible only 
for the borderline values of the parameters, 

H ~ e 1/4 ~ a 1/3 ~ KT 3 eV. (46) 

However, as we discussed at the end of Section IV.A, these values seem to be disfavoured 
by observations. 

The brane nucleation rate at present is given by Eq.(^) with B and A from Eqs.(|9|),(|Tl|). 
In order to have no brane nucleations in the observable universe in a Hubble time, we have 
to require that 

r*o - 1. ( 47 ) 

where to is the present cosmic time. For the parameter values fl4"BD , A ~ (10 -3 eV) 4 and 
Eq.(g7|) gives exp(-B^) < KT 116 , or 

B (flat) > m ( 4g ) 

This is only marginally consistent with the bound (|25|) 

The brane nucleation rate during inflation is determined by the smaller of the two bounce 
actions (|9l), (|T0l) . Eq.([5J) tells us that in models based on M-theory brane nucleation can be 
suppressed only if the expansion rate during inflation is H ~ 10~ 3 eV. Let M ~ 30 be the 
number of e-foldings from the time when the comoving region corresponding to the presently 
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observable universe crossed the horizon to the end of inflation. Then the size of this region 
at the end of inflation is H^e^. In order to have no brane nucleations in this region during 
this whole period, we have to require 



YH^e m < 1. (49) 



For the parameter values (fig) this gives B £ 90, again marginally consistent with (pq) 

We thus see that M-theory based four-form models could in principle provide a solution 
to the cosmological constant problems, but only if inflation is at a TeV scale and a and e 
are in the tight corner of the parameter space (|46|) . With such values of the parameters, the 
condition (fig 1 ) can be (marginally) satisfied. However, from (|TJ] we then find a significant bias 



towards A-lowering nucleation events, which would shift the a priori distribution (^) towards 
lower A. This would result in a prediction near the lower anthropic bound p\ ~ — Pmo- The 
bias towards A-lowering events might be compensated to some extent by the differential 
expansion rate, which adds relative volume to regions with high A. However, both effects are 
exponential, and unless there is a conspiracy in the parameters of the model, the differential 
expansion is likely to be either insignificant or dominant. In the latter case, the a priori 
distribution would be biased towards large A, and it would be likely to predict a cosmological 
constant much larger than observed. In summary, it seems difficult to obtain a flat a priori 
distribution even in the range (46). Of course, the possibility cannot be excluded with our 
order of magnitude estimates, and there may still be a small viable region of parameter space 
in this borderline range. We note also that for models unrelated to M-theory the allowed 
parameter space is much larger. 

Suppose now that the visible universe contains more than one bubble. This would 
generally result in microwave background anisotropies of amplitude 5T/T £ e/ Pmo, so to 
avoid conflict with observations we have to require 

e & 10- W (50) 

This takes us far from the borderline values (fig 1 ), and thus the multiple bubble scenario 
cannot be realized in M-theory based models. For non-M-theory models, a suitable set of 
parameters can be easily found by choosing a and H sufficiently large, while keeping e under 
the bound ©. 

The multiple bubble scenario is feasible only if branes have negligible interaction with 
ordinary matter. Otherwise we would see fireworks along the bubble boundaries, where the 
branes hit the stars and where they hit one another. However, the gravitational impact 
of the branes cannot be avoided. An observer outside an expanding spherical bubble does 
not experience any gravitational force until he is hit by the brane. While the brane passes 
through the observer, the part of his body inside the brane will experience an acceleration 
a = GM/R 2 relative to the part of the body still outside the brane. Here, M = (4n/3)eR 3 
and R is the bubble radius at the moment of impact. With R ~ t and e satisfying (0), we 
have a ~ Get ~ (e/ Pmo)^o 1 ~ 10~ 12 cm/ s 2 . The relative speed developed during the passage 
time At ~ 10 -8 s is Av ;$ 10~ 20 cm/s, and the corresponding displacement is much smaller 
than the inter-atomic distance. For a brane passing through a Sun-like star, At ~ 10 s and 
the displacement is still smaller than the atomic scale. Thus, if a brane is to sweep through 
the Solar system, its only effect would be to set up imperceptible vibrations in the objects 
it leaves behind. 
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What happens if B < 1, so that brane nucleation is completely unsuppressed? The main 
danger here is that the vacuum energy will decay so fast that it will drop significantly in 
less than a Hubble time. This can be countered by choosing e so small that the change in 
Pa is negligible even after nucleation of a large number of bubbles. This case, however, is 
almost indistinguishable from that of a scalar field with a very flat potential, which will be 
discussed in Section VI. 



C. No empty universe problem 

Here we shall comment on the so called empty universe problem which was encountered 
in all earlier work on discrete A models [31,19 J25|j5|,|^,ptp6l . The scenario these authors had 



in mind is that the universe starts with a large cosmological constant and relaxes, within 
the available cosmic time, to a metastable vacuum with an observationally acceptable value 
of A. The problem is that, in order to make the present vacuum sufficiently stable, the 
brane nucleation has to be strongly suppressed. One then finds that the time it takes the 
universe to evolve to the low-energy vacuum is so large that, by the time when the process 
is complete, any matter that the universe initially had gets diluted to an extremely low 
density. So one ends up with an empty universe dominated by the cosmological constant. 

A number of solutions to this problem have been proposed. FMSW suggested that the 
nucleation rate of multiple coincident branes may be enhanced due to the increased density 
of states. They argued that this would lead to a rapid descent of the vacuum energy towards 
lower values. To ensure the long lifetime of the present vacuum, they argued that this rate 
enhancement may not apply to the vacuum with the lowest positive value of p\. Bousso 
and Polchinski ||, who considered brane nucleation with large jumps in pa, suggested that 
the penultimate vacuum could have a high energy density. The inflaton field would then be 
excited to high values of its potential by quantum fluctuations. When the ultimate brane 
nucleates, the inflaton rolls down the potential thermalizing its energy and providing a high 
density of matter. Alternatively, they suggested that the nucleation of the ultimate bubble, 
which in their model is accompanied by a large change in the four-form field F, can be 
accompanied by a large modification of the inflaton potential. As a result the inflaton will 
be displaced from the minimum of the potential, even if it was at the minimum prior to the 
bubble nucleation. 

In our view, the empty universe problem is not a real problem, and the attempts to 
solve it seem therefore unnecessary. The problem disappears when the eternal nature of 
inflation is taken into account. As the inflaton fluctuates back and forth in the quantum 
diffusion regime, branes are constantly being nucleated and all possible values of p\ are 
reached. The slow rate of brane nucleation is not a problem, since an unlimited amount of 
time is available. Thermalization of the inflaton energy occurs at different times in different 
parts of the universe, and each region inherits the local value of p\. Each possible value 
is represented in the thermalized regions of the universe. We are interested only in those 
regions where p\ is in the anthropic range (H), because that's where all the galaxies are. 
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V. MULTIPLE BRANE NUCLEATION 



Up till now we assumed that brane nucleations change the four-form field F by a single 
unit, Eq. (|20|) . However, nucleation of multiple coincident branes is also possible. For k 
coincident branes there is a U(k) super Yang-Mills (SYM) living on the world-volume. In 
FMSW ||] it was argued that the nucleation of coincident branes would be enhanced by a 
large degeneracy factor 

D = e 5 , 

where S is the "entropy" of the SYM fields. For 2-branes arising from the wrapping of a 
4.D-brane on a degenerating 2-cycle, FMSW estimated this entropy as 

S ~ k p R 2 T 2 . (51) 

Here k 13 counts the effective number of degrees of freedom which live on the brane. There are 
theoretical uncertainties in the exponent (3, but FMSW suggest that it should be between 
2 and 3/2. R is the radius of the Coleman-De Luccia instanton, which coincides with the 
size of the "bubble" at the time of nucleation, and T is some effective temperature. FMSW 
considered two different candidates for the effective temperature. One of them was the 
effective ambient de Sitter temperature [^2j T before brane nucleation, and the other was 
the geometric mean of T and the effective temperature Tj of the new de Sitter space inside 
the nucleated brane T ~ (TjTo) 1 ^ 2 . 

It is easy to understand, however, that the relevant effective temperature corresponding 
to the Coleman-De Luccia (CdL) instanton is in fact none of the above, but simply the 
effective de Sitter temperature of the 2+1 dimensional world- volume of the brane 

T = — . (52) 

This is the temperature experienced by the degrees of freedom living in the wall (and it is 
in fact higher than T and Tj). The prefactor D is a determinant arising from Gaussian 
integration of perturbations around the instanton solution, including all light degrees of 



freedom. Such determinants where discussed in some detail in 231. A scalar field of mass 



m living on the 2+1 dimensional world- volume gives a contribution 

D s = e« 2 (53) 

where ((z) is the Zeta function of the scalar fluctuation operator on the 3-sphere. Its 
derivative at the origin is given by 



2 r^y 

('(0) = 2(' R (-2) - y 2 ln(sin iry) + — / x ln(sin x)dx, (54) 

n 2 Jo 

where y 2 = 1 — m 2 R 2 and C,r is the usual Riemann Zeta function (this expression is valid 
for light fields, with mR < 1). 

For instance, the contribution of a conformally coupled scalar field can be obtained by 
taking m 2 = (3/4)R~ 2 , which gives 

C'(0) = 2&(-2) - ^ln2 + -^C*(3) « -0.1276 

4 S7T Z 
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Hence, the effective degeneracy factor contributed by a conformal scalar field is given by 



D conf . « e- 0638 « 0.91 < 1. (55) 

The first thing to note is that this factor is not an enhancement, but a suppression. Hence, 
the determinant cannot simply be thought of as the exponential of an entropy. 

In fact, the CdL instanton is not a thermal instanton, but a zero temperatue instanton. 
Thermal instantons for brane nucleation are static and have the topology S 2 x S 1 (rather than 
S* 3 ), where the S 2 is the 2D-brane at fixed time and the S l is the periodic Euclidean time. 
Thermal instantons do in fact exist also in de Sitter space, but they have not received too 
much attention because their Euclidean action is always larger than that of the maximally 
symmetric CdL. For thermal instantons (in flat or in de Sitter space) the determinantal 
prefactor is given by D = e~ AF/T , where AF = AE - TS is the free energy contribution of 
light degrees of freedom on the brane. Such prefactors have been considered in ||33|| . The free 
energy consists of the vacuum energy AE (or Casimir energy on the two-sphere) minus the 
product of the temperature times the entropy. While the entropy is always positive, the sign 
of the Casimir contribution is notoriously dependent on the type of field. In fact, for thermal 
instantons in de Sitter the temperature is always smaller or equal to the inverse of the size 
of the bubble, and hence the sign of the free energy contribution can easily be dominated by 
the Casimir contribution, which can have either sign. Although as mentioned above the CdL 
instanton is not thermal, this consideration may help clarify why the prefactor D need not 
represent an enhancement. Depending on the field content it may represent a suppression. 

Another thing to note about ([55]) is that it is independent of R. Roughly speaking, this 
is consistent with the fact that the effective temperature is T cx Br 1 . In general, however, 
the degeneracy factor will depend on R and on the mass of the particle. For light minimally 
coupled scalars, equations ( pB|) and (^) give 

e Ca'(-2) 

D '~^R < mfl<<1 >- 

There can be a strong enhancement in the nucleation rate if there are very light massless 
scalar fields. In the limit m — > the factor goes to infinity. This is because a massless scalar 
has a normalizable zero mode on the sphere, corresponding to the symmetry <p — > <p + const. 
In this case, the zero mode must be treated as a collective coordinate. The nucleation rate 
is proportional to the range 8(j> of the field 0, because the bubbles can be nucleated with 
any average value of the scalar field with equal probability |23[ 



5D s (m 2 = 0) = lim [m J D s (m)](vr J R 3 ) 1 / 2 (50 = e^'^R 1 ^. 

As we shall discuss in some examples below, some scalars are likely to pick up masses of 
order of the intrinsic curvature of the 2+1 sphere, and for these D s is also independent of 
the radius. 

Let us briefly consider the field content on the brane. For k coincident 4D-branes in 
ten dimensions, the effective theory is U(k) SYM. This consists of a U(k) gauge field plus 
5(k 2 — 1) scalar degrees of freedom in the adjoint representation of SU(k) plus 5 scalar 
singlets plus the corresponding fermionic degrees of freedom. 

If the branes are flat (as in the case when there is no external four- form field), then the 
theory is supersymmetric and all scalar degrees of freedom are massless. For the case of 
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a single brane, the five scalars represent the goldstone modes of the broken translational 
symmetry. That is, they correspond to transverse displacements of the brane. For the case 
of two branes, there are 10 such displacements. Five of them correspond to simultaneous 
motion of both branes. These are the singlets under SU(2). The rest are in the adjoint 
representation, and if they acquire an expectation value they give mass to two of the four 
gauge bosons. For instance, when the two branes move apart, one of the adjoint scalars 
acquires an expectation value and two of the gauge bosons get a mass, breaking the symmetry 
27(2) -»• 17(1) x 17(1). 

The case of interest to us is not a flat brane, but a 4/J-brane wrapped on a degenerating 
two cycle. The world-wolume of the resulting 2/J-brane in 4 Euclidean dimensions is not 
flat either, but forms a 3-sphere of radius R. In this situation, we do not expect the theory 
to be super symmetric. This is just as well, since in order for the instanton to make any 
sense, some of the adjoint scalars must pick up masses at one loop. Otherwise the instanton 
would have too many zero modes and too many negative modes. To simplify, let us consider 
the 3-sphere in 4 non-compact dimensions. For the case of a single brane, the transverse 
displacements correspond to a scalar field of mass m 2 = —3R~ 2 |[34|| . This scalar field has 
a single negative mode, which is the constant I = mode. A negative mode is precisely 
what is needed for an instanton to contribute to the imaginary part of the vacuum energy, 
and hence to contribute to false vacuum decay Also, there are four normalizable zero 



modes, which are the spherical harmonics with 1 = 1. These correspond to the four space- 
time translational modes of the instanton, which have to be treated as collective coordinates. 



This Goldstone field gives a determinantal prefactor of the form |23 



D = ^e^(" 2 »fi, 

where Q = VT is the spacetime volume. The prefactor A in the nucleation rate (|Tl~D is 
obtained after dividing by Q. 

If there are 2 coincident branes, then in principle there would be two such fields 0i and 
2 corresponding to the independent transverse displacements of the brane. However, only 
the combination + = (0i + 02) will correspond to a singlet under 577(2). The orthogonal 
combination 0_ = (0i — 02) will be in the adjoint. As mentioned above, if the instanton 
with two coincident branes is to make any sense, this combination must aquire a positive 
mass at one loop so that there is a single negative mode, not two, and four normalizable 
zero modes in total. In other words, in order for the instanton to make sense, the branes 
must attract each other. If they repelled each other or if they did not "interact", then the 
two brane configuration would in fact be an accidental superposition of two independent 
bubbles in the "dilute gas" of instantons. The mass of 0_ can be estimated as follows. 
The mass of the gauge field is given by m^(0-) ~ M P f(d), where M P is the Planck 
mass, d is the distance between branes and f(d) ~ d for d ^> Mp 1 . This is because this 
vector corresponds to fundamental strings stretching from one brane to the other. For 
smaller distances, d ^ Mp 1 , we may expect a milder behaviour for the mass, which we may 
heuristically parametrize as a power f(d) ~ d(Mpd) n , with n > 0. The canonical field is 
related to the distance through 0_ ~ da 1 / 2 . Hence, 

ml-Mpfia- 1 ^). 
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On a flat brane, the effective potential induced by a gauge field of mass tua is proportional 
to m\. However, it can be shown that on a sphere there is also a term of order m\R~ l which 
will in fact dominate at small m^. When these terms are added to the tree level potential 
— Si? -2 */) 2 , the scalar aquires a very tiny expectation value < </>_ >~ (a/ RMp) 1 / 2n (T 1 / 2 Mp 1 . 
In the broken phase, </>_ has a positive mass squared of order m\_ ~ nR~ 2 . The gauge 
bosons will in turn aquire imperceptibly tiny masses m\ <C R~ 2 . 

To summarize, some of the scalars may get very large masses from the wrapping on a 
degenerating cycle. These will decouple, and presumably will not contribute to the degener- 
acy factor. Others, corresponding to the relative position of the 2-branes in 3+1 dimensions 
will get masses of order -R -1 , and hence will contribute degeneracy factors of order one, just 
like the conformally coupled field discussed above. For the massless or nearly massles gauge 
fields the contribution to the degeneracy factor will be of order one (the vectors have no 
zero modes on the sphere, so unlike the case of a scalar, a tiny mass will not cause a large 
degeneracy factor). Similar considerations could be applied to fermions. Thus we expect 
the total prefactor to be of order 

D total ~ tta 2 R 2 e akP (56) 

where a is a numerical factor and k@ is, as in Eq. (|5"ID , an estimate of the effective number 
of degrees of freedom. For a flat brane at weak coupling, (3 = 2, but as argued by FMSW 
it could be lower for the wrapped brane. Unfortunately, without going into a very detailed 
analysis (which is outside the scope of this paper) we are unable to determine the sign of 
the constant a. However, as argued above, this value seems to be rather insensitive to the 
value of R or to the value of the ambient de Sitter temperature. 

For a < the nucleation of multiple branes is suppressed and we are back to the situation 
described in Subsection IV. B. For a > a disaster may occur because transitions into deep 
Anti-de Sitter space through multiple brane nucleation seem to be unsuppressed due to a 
large degeneracy factor. FMSW suggested that the disaster could be averted by an anthropic 
argument. If the step e in the vacuum energy is of the order of 3 in units of Pmo, allowing the 
values ...,-2,1,4,... then the stringent anthropic bound for a negative cosmological constant 
Pa ~ _ Pa/o would tell us that the vacuum energy is in fact the lowest allowed value in 
the list, that is pa ~ pmo (note that this argument requires a certain adjustment of the 
step e in order to explain the observed value). However, there is another problem which 
is how to explain the stability of this vacuum once it has been reached. In the FMSW 
scenario, the stability was attributted to the fact that in the vacuum with pa ~ Pmo the 
effective temperature of the brane would be so low that the degeneracy factor is switched 
off. However, as we have seen, the degeneracy factor is quite independent on the ambient 
de Sitter temperatures and hence it does not seem to switch off. The same mechanism that 
would enhance coincident brane nucleation from a high energy vacuum, would cause the 
disastrous decay of "our" vacuum. 

Finally, we note that even if a mechanism could be devised to switch off the degeneracy 
factor, so that the present vacuum is stable, the time coincidence t G ~ t\ would be left 
unexplained by this approach (just as in the non-anthropic models discussed in Section 
VIII). Also, the unsuppressed nucleation of coincident branes seems to preclude eternal 
inflation, and even if one may intuitively argue that the lowest anthropic value is the most 
probable, the actual probability distribution for positive pa seems hard to estimate. 
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VI. SCALAR FIELD WITH A VERY FLAT POTENTIAL 



In this class of models, what we perceive as a cosmological constant is in fact a potential 
V(4>) of a scalar field <f)(x). The potential does not have a succession of minima as in Abbott's 
washboard model, but its slope is assumed to be so small that the evolution of is slow on 
the cosmological time scale. This is achieved if the slow roll conditions 

MpV" < V < p M0 , (57) 

M P V < V < pa/0, (58) 

are satisfied up to the present time (here it is assumed that any "true" cosmological constant 
is also included in the potential V .) These conditions ensure that the field is overdamped by 
the Hubble expansion, and that the kinetic energy is negligible compared with the potential 
energy (so that the equation of state is basically that of a cosmological constant term.) The 
field is also assumed to have negligible couplings to all fields other than gravity. 

Let us now suppose, as in the previous sections, that there was a period of inflation 
driven by a different scalar field x- During inflation, massless scalar fields are randomized 
by quantum fluctuations, which cause their root mean squared value to increase with time 
as A0 ~ H(Ht) 1 / 2 , where H is the inflationary expansion rate. If we consider a field of 
mass m, this effect competes with the classical drift down to the bottom of the potential, 
and after some time of order t ~ HmT 2 a stationary distribution with root mean squared 
A0 ~ H 2 m~ l is established. This can be interpreted in terms of the Gibbons-Hawking 
temperature T ~ H of de Sitter space as the condition V{<p) ~ m 2 2 ~ T 4 . In this example, 
all field values |0| <C H 2 jm would be almost equally probable after the end of inflation. This 
discussion, however, assumes that inflation proceeds at almost the same rate for all field 
values in the range considered. That is, the differential expansion rate 6H ~ V (<p)(H Mp)^ 1 
is ignored. 

The case of interest to us is slightly more general because the potential need not be 
quadratic, and also because we are not necessarily interested in field values near = 0. 
Rather, we are interested in field values for which the energy density is in the anthropically 
allowed range 

- PA/o £ V{4>) < lOOpA/o- (59) 

The differential expansion rate 5H ~ V (<f))(H Mp)~ l will be negligible if the time t ~ 
(A(j)) 2 nth H~ 3 that it takes for the field to fluctuate accross the anthropic range of corre- 
sponding to (|59| ) is smaller than (SH)" 1 for the same range. This requires 

H 2 

(A(j)) anth < m M P . (60) 

IOPa/o 

If this condition is not satisfied, then the a priori probability for the field values with a higher 
V(0) would be exponentially enhanced with respect to the field values at the lower anthropic 
end. This would result in a prediction for the effective cosmological constant which would 
be too high compared with observations. Therefore, in what follows, we shall demand that 
our potential satisfies Eqs. Q57|) , (]58[) and (|60|). 
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A. Solving both cosmological constant problems. 



Consider a potential of the form M, 

V{4>) = Pbare + \p 2 <t>\ (61) 

where pbare represents the "true" cosmological constant. If pbare and p 2 have opposite signs, 
then the effective vacuum energy will be very small when 

|0| « \2 Ph are\ ll2 \pl l - (62) 

The anthropic range is given by (A(p) anth ~ 100pMo|/^ 2 Pbare| _1/ ' 2 - Then, conditions (0), ( |58|) 
and ( |S0D imply 

1/2 

10 H 2 \p ba re\^M P ^ ™ ^ \ Pha re\^M P ^ 

From the cosmic microwave background temperature fluctuations we know that 
H £ 10 _5 Mp. This leaves a wide range of possibilities for the value of the mass param- 
eter, 



Ml 

\Pba 



/^(lO-^-lO- 120 )^^^, (64) 



spanning some 47 orders of magnitude. Provided that p is in this range, the a priori 
probability distribution V*(<fi) for will be flat. The probability distribution for the effective 
cosmological constant p^ = V((f>) is given by 

and it will also be very flat, since V is almost constant in the anthropic range. As men- 
tioned in Section II, a flat a priori distribution for the effective cosmological constant in the 
anthropic range entails an automatic explanation for the two cosmological constant puzzles 

00. 



B. A small mass from instantons? 

The challenge in the scenario presented above is to explain the small mass parameter ( |64|) . 
In Ref. 0] we suggested that this can be achieved through instanton effects. For instance, 
4> could be a pseudo-Goldstone, the phase of a scalar field $ = rje 1 ^^ which spontaneously 
breaks a global U(l) symmetry. Since global charge can be swallowed by wormholes, a small 



mass term for the field <p will be generated through gravitational instantons [35,36 



Another possibility is that the phase may have an "axion" coupling of the form 

—(f)FF, (65) 
V 
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where F is the field strength of a "hidden" gauge sector with gauge coupling constant a s . 
The coupling (|65|) will give a small mass to the pseudoscalar through instanton effects. 

We should mention, however, that there may be certain limitations in implementing this 
idea in the present context. Consider an instanton-induced potential of the form 

V{4>) = Pbare + A 4 COs(0/^) . 

In order to solve the cosmological constant problem we need 

A Pbare- 

Combining this with the slow roll conditions ( jSffi ) and (|58p we find 

ji r Pbare 

r] > M P . 

Pa/o 

Thus, the expectation value rj must be truly huge compared with the Planck scale. 

In usual axion models, the effective vertex (|65|) can be obtained in the following way. The 
scalar field $ has Yukawa interactions of the form h&^ty with an "exotic" fermion \l/ (here 
h is the Yukawa coupling constant). The fermion in turn interacts with the non abelian 
gauge fields, and the coupling flBSp is generated at one loop. The mass of the fermions in 
the broken phase is given by ~ hrj. In our case, this mass is extremely large (unless h 
is extremely small), and so we can hardly trust the field theory model for generating flB"5T). 

Perhaps more worrisome is the effect of wormholes. For small symmetry breaking scale 
rj ,iS Mp the scale A 4 in the instanton potential is of order Mp e~ s , where S ~ Mp/rj is the 
wormhole action |35| , |36| . The radius of the wormhole is given by R ~ [Mpr])- 1 / 2 . This radius 



approaches the Planck scale as 77 approaches M P , and the process becomes unsuppressed. 
The instanton calculation becomes unreliable for higher values of rj, but it is not clear what 
would prevent nonperturbative gravitational effects from completely destroying the global 
symmetry. 

Therefore, as mentioned above, the generation of a small mass through instantons may 
not have a straightforward implementation in the present context. Clearly, this issue deserves 
further investigation (see e.g. |36[|). 



C. A very flat potential from field renormalization 

Consider a potential of the form |1| 

V(0) = pw e + M 4 /(A0), (66) 

where M is a reasonable mass parameter and / is a function of order one with no large or 
small parameters. If the parameter A in the argument of / is chosen to be very small, then 
V(4>) will be very flat. In particular, the mass term of the field which has two powers of A 
will be very small. Weinberg suggested [0 that perhaps the smallness of A can be attributed 
to a large running of the field renormalization from some fiducial short distance scale p 
to the large scales in which the cosmological constant is relevant, p —>■ 0. 

More generally, as noted in Ref . , the effective Lagrangian for a scalar field ip at large 
distances will include nonminimal kinetic terms, 
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C = F\ij){d^f-V^) + 



(67) 



Here, F plays the role of a field renormalization, which in fact may depend on if), and we 
have omitted terms with more derivatives of if). If F is very large, then the field redefinition 

dip = Fdip 

will result in a very flat effective potential for <f>. 

Take for instance F = e^' Mp and V = Pbare + (m 2 /2)ip 2 , where m is a not too large mass 
parameter (see below). After the change of variables we obtain 

C = (d^) 2 - Pbare - ~m 2 M 2 P [HMMp)] 2 + ... 

The effective potential is now very flat at large <f>. The slow roll conditions (|57| ) and ( |58| ) 
are satisfied for 

, . , , . m 2 m 

The antropic range m 2 if) 2 ~ |pf, are | will satisfy this condition provided that m <C 
IPbare^^Mp 1 . Finally, the condition flBTf ) is easily satisfied by chosing a sufficiently high 
Hubble rate during inflation H 2 3> lO 3 pQ0(0/0 m j n ). 

Thus, starting from a Lagrangian fl6?D with fairly simple functions F and V we have 
been able to satisfy all necessary conditions to solve both cosmological constant problems. 
Of course, one may wonder why F should have exponential behaviour when V is only poly- 
nomial, and it would be good to find a well motivated physical setup where this Lagrangian 
emerges in a natural way. 



VII. A SLOWLY VARYING FOUR-FORM? 

In theories with extra dimensions, the four form field strength is dynamical above the 
compactification energy scale. Donoghue suggested that in the early universe the four- 
form might take a continuum of different values in different parts of the universe, and that 
it might get frozen to these values as the universe cooled down below the compactification 
scale. However, it is easy to show that the effective cosmological constant can vary from 
place to place only if the size of the internal space is also variable. As a result, the effect of 
the four form is more properly described as a contribution to the effective potential for the 
radius modulus of the extra space. 

Let the higher dimensional manifold be the product of a four-dimensional spacetime M. 
and an internal space S, 

ds 2 = g iliV [x)dx il dx v + a^{x, y)dy l dy 3 ' , (68) 

where \i, v — 0, 3 are the four-dimensional indices and i,j = 1, n are the internal space 
indices. The field strength takes the form, F = f(x,y)u)4 + where uj± = ^fg (A^dx^) is 
the four-dimensional volume form (g = — detg^) and the ellipsis denote terms with at least 
one internal index (these do not behave as a four- form upon dimensional reduction). The 
equations of motion reduce to 
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dG = 0, (69) 

where G =* F = f(x,y)u n + ... is the dual of the field strength and u> n = yfa (Aidy 1 ) is 
the volume form on the internal space (a = deta - ^). Consider two different points x% and 
x 2 on the 4D manifold Ai, and a curve 712 joining them. Applying Stokes theorem to the 
"cylinder" 712 x S (where S is the internal space), and using the equations of motion (|69|), 
we immediately find that 



G - / G = / dG = 0, (70) 

S;xi J S]X2 JSXJ12 

for arbitrary x\ and x-i- Assume that the internal metric factorizes as 

<ry = e^a^y). 

Then, Eq. ( |70f ) implies 

fix) = f e™H*\ 

where f(x) is the average of f(x,y) over the internal space and fo is a constant. Kaluza- 
Klein modes average to zero on the internal space and do not contribute to /. However such 
modes are massive in the reduced theory and do not behave as an effective cosmological 
constant. It follows that the contribution of the four form to the effective cosmological 
constant is 

\F 2 = ^J s F^F^°^ d n y = 1f 2 c -*« (71) 

where Fq = /q J s \J~d d n y = const. 

In the dimensionally reduced theory, ip(x) is a four dimensional scalar field, and (|7T|) is 
just a contribution to its effective potential V(ip). At the classical level, there are two other 
such contributions, due to a bare higher dimensional cosmological constant and due 

to the curvature of the internal manifold. Following Refs. [37], it is easy to show that in 



terms of the Einstein frame metric g^ v = e'^g^, the effective action takes the form 



R K — — '- g^d^d^ - V(if>) 



with 



Vty) = A e~ n ^ + ^ e" 3 "^ - j e ~^ 2+n ^ . 



Here Mj> = M 2+n Vo and A = A^^Vq, where is the higher dimensional Planck mass, 
— Is V&d n y = const., and K is defined by = (K/n)<Jij. In general, the potential may 
have a minimum but this need not be near V = 0. For n > 1 one can adjust the parameters 
A and F so that the minimum is at the right height to fit observations , but this would 
be the usual fine-tuning (for n = 1 the minimum must have negative effective cosmological 
constant, so this tuning is not possible). We may also consider the possibility that the field 
ip is away from the minimum, but slowly rolling so that the effective potential V(ip) plays 
the role of an effective cosmological constant, as described in Section VI. The problem is that 
if the slow roll conditions are met, then ip would have a negligible mass and would mediate 
long range interactions of gravitational strength, which are ruled out by observations. 
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VIII. NON-ANTHROPIC APPROACHES 



Here we comment on some attempts to solve the cosmic coincidence problem without 
resorting to the anthropic principle. In a recent paper [|7| Arkani-Hamed et. al. suggested an 
explanation to the approximate coincidence of several cosmological timescales: the time of 
matter-radiation equality t eq , the time of A-domination t\, and the time of galaxy formation 
to- They assume that the Planck scale Mp and the electroweak scale M w are the only 
relevant scales and argue that the temperature at matter-radiation equality and the vacuum 
energy should then be given by 

T eq ~ Ml/M P , (72) 



p A ~ {MHM P f. (73) 



It follows immediately from (|72"|) , (|73"1) that t eq ~ t\. This coincidence should of course be 
understood in a very rough sense, since the actual values of t eq and £a in our universe differ 
by a few orders of magnitude. Now, assuming the density fluctuation amplitude 

Q = 5p/p ~ 10~ 5 , (74) 

and using a more accurate value for t eq /t\, the authors show from (|73|) that the epoch of 
galaxy formation is at 

to ~ t A . (75) 
In our view, a relation like (|73|) may account for the smallness of A and may even explain 



its observed value. However, the cosmic coincidence (|75j) would remain unexplained. The 
time of A-domination is determined by the value of A, while the epoch of galaxy formation 
is determined by the amplitude of density fluctuations Q. Even if we explain the value of A, 
we still have to explain why the value of Q is such that t G ~ t\. Moreover, the accuracy of a 
few orders of magnitude is not sufficient to explain the cosmic time coincidence: observations 
indicate that the coincidence flT5| ) is accurate within one order of magnitude. 

Another non-anthropic approach to solving the cosmic coincidence problem involves k- 
essence, a scalar field with a non-trivial kinetic term ||. fc-essence has a positive effective 
pressure during the radiation era and starts acting as an effective cosmological constant 
with the onset of matter domination. With a suitable choice of parameters it dominates 
the universe at t\ ~ to- However, one could also choose parameters to obtain t\ <C to or 
t\ ^> to- This model can explain why A-domination occurs at t > t eq , but it cannot account 
for the coincidence ([75]). 



IX. MODELS WITH VARIABLE Q 

Several authors have recently expressed the view that the anthropic principle can perhaps 
be applied to the cosmological constant problem - but to nothing else! [§^,0 For instance, 
Weinberg has remarked [|18|] that we cannot explain the masses and charges of the elementary 
particles by assuming that they depend on the expectation values of scalar fields with very 
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flat potentials. These light fields would couple to the elementary particles, and would have 
been observed in collisions and decays. 

While this remark may be true, we can still apply the anthropic principle to variables 
which determine the large scale properties of the Universe, and which generically fall into 
the category of "initial conditions". Examples of these are the amplitude of primordial 



fluctuations Q |4ll Jffijlql , the density parameter f2 [^0| , p9f| , or even the baryon asymmetry. 
In the inflationary context, these parameters depend on the path that the inflaton field takes 
in going from the diffusion regime to thermalization. The inflaton potential represented in 
Fig. 2 is one-dimensional, and there is a single path from the top of the potential to the 
local minimum. However, more generally, the inflaton has several components, and there 
may be a continuum of paths from the diffusion region to a given minimum. Even if the 
low energy particle physics Lagrangian is the same in all thermalized regions, and even if 
there are no exotic light degrees of freedom after thermalization, these regions may start 
with different initial conditions which will be more or less favourable to galaxy formation. 

Consider for instance |39] a two component scalar field x = Xi + ^X2 — \x\ el@ i with 
potential V(x) = (91X1 + 92X2)/^- This potential produces inflation for |%| ^ Mp. However, 
the amplitude of density perturbations Q depends on the direction G of approach to the 
minimum, Q ~ m(Q)N(\x\)Mp l . Herem 2 (6) = gi cos 2 6+^2 sin 2 6 and iV ~ \x? M p 2 ~ 60 
is the number of e-foldings from the time the present Hubble scale first crossed the horizon 
until the end of inflation. The minimum at x = will be reached from different directions 
in different thermalized regions, and therefore these regions will have a different value of Q 
as an initial condition. This example illustrates that Q can easily be made into a random 
variable. In general, its a priori distribution V*(Q) (i.e. its volume distribution at the time 
of thermalization) will not necessarily be flat in the anthropically allowed range. For any 
given model, this distribution can be calculated using the numerical methods of Ref. |2U . 
To proceed, however, we shall heuristically parametrize it as 

dV^Q) ~Q~ a d\nQ, (76) 

where a is a constant (this may not necessarily be a good estimate for the particlular model 
presented above, but we shall use it anyway for the sake of argument.) 

We may now take a point of view which is "complementary" to the one used in the 
preceeding sections. Let us assume that the cosmological constant is truly a constant of 
order M®Mp 6 determined from the fundamental theory (as assumed e.g. in 0), and that 
Q is a random variable with prior distribution (|76|). If a > 0, then low values of Q will be 
favoured a priori. However, if Q is too low, galaxies will not have time to form before the 
time t\ when the cosmological constant starts dominating. With this, we would basically 
explain why Q ~ 1CT 5 as well as the time coincidence t G ~ t\. These arguments can be 
made more quantitative in the following way. The probability distribution for a galaxy to 
form at time t G is given by 

dV(t G ) oc V*(Q) MtG ' tA ' Q) d\nQ dt G . (77) 
at G 

Here, v{t G , t\, Q) is the fraction of matter that clusters up to the time t G in a universe where 
the density contrast at the time of recombination is Q and where the cosmological constant 
is such that it will start dominating at the time t^. This fraction can be easily estimated 
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by using the Press-Schechter approximation. In Ref. [Tj| it was shown [see Eq. (27) of that 
reference] that after integrating over Q the probability distribution for to is given by 



dF° 

dV(t G ) oc — — dx 



dx 



(78) 



where 



and 



F(x) 



1 + x 



X 



1/2 



dw 



o w^il+w) 3 / 2 ' 



X 



sinh 2 (t G /t/ 



(Following [16], we are using the convention that £a is the time at which £7a = ^Msinh 1, 
where Qi are the fractional densities of cosmological constant and non-relativistic matter 



respectively.) The distribution (|78|) is plotted in Fig. 3 for different values of a, and we see 
that for moderate values of a it presents a rather prominent peak at t G ~ t\, as expected 
from the general arguments above. 




FIG. 3. Probability distribution for tc/tA, the time of galaxy formation as compared to the time 
when the cosmological constant starts dominating. Here, A is taken to be a fundamental constant 
but the density contrast Q is treated as a random variable with a priori volume distribution oc Q~ a 
at the time of thermalization. The plot is shown for a = 1.5, 3 and 5. The distributions present 
rather well defined peaks at to ~ £a- 

Finally, one may take the view that both A and Q are random variables. This possibility 



was considered in [16], where it was shown that a decreasing a priori distribution for Q 



pushes the cosmological constant to small values, so that both tc and t\ tend to be very 
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large. In this case, a new time scale comes into play. This is the so-called cooling boundary 



td, ||41|| . For times t > gravitationally collapsing clouds of galactic mass cannot fragment 
into stars because they are too cold to reach the usual "cooling" line emission thresholds, 
and they stay as pressure supported configurations for a very long time. Thus, usual galaxy 
formation is suppressed after t ~ t c b ~ 3-10 10 Fr. This time is determined from microphysical 
parameters such as the fine structure constant, the proton mass and the fraction of baryonic 
matter. Since the time of galaxy formation cannot be arbitrarily large, in the situation 
where both Q and A are random variables we expect t G ~ ~ t cb (see [|16| for details). 



There are many uncertainties associated with the calculation of t^. Perhaps after some of 
these uncertainties are removed, we may actually find that ^> t G ~ ^a- This hypothetical 
situation would suggest that one of the timescales t G and t\ is not a random variable, or that 
if both of them are, then their a priori distribution must have a rather peculiar behaviour. 
This in turn would give us information on the theories of initial conditions giving rise to 
these a priori distributions. 

These examples seem to suggest that the applicability of anthropic reasoning, once it is 
accepted, may easily go beyond the issue of the cosmological constant problem. 



X. CONCLUSIONS 

The anthropic principle has a bad reputation. It is often regarded as a handwaving argu- 
ment relying on poorly understood phenomena like intelligent life and having no predictive 
power. Although this criticism is not entirely ungrounded, there is a class of cosmological 
models where the use of anthropic principle is not only justified but may in fact be in- 
evitable, and where it can be used to make quantitative predictions. These are the models 
in which some cosmological parameters, or physical "constants", take different values in 
different parts of the universe. In such models, one cannot predict the precise values of the 
parameters that we are going to observe. One can only hope to calculate the corresponding 
probability distributions. The criteria for justifying (and compelling) the use of anthropic 
principle are that the model should provide (i) a mechanism for variation of the parame- 
ters and (ii) a way of calculating the probability distributions. Once the probabilities are 
calculated, one can predict that the parameters are going to be observed within a certain 
range of values, say, at a 95% confidence level. This seems to be as quantitative as one can 
possibly get in this class of models. 

From a practical point of view, parameters that we can hope to determine anthropically 
should satisfy the condition that they do not affect life processes, and preferably should also 
not affect the poorly understood astrophysical processes such as star formation |^ , |TE| , |2"6" 



So for example the gravitational constant may be hard to determine anthropically with our 
present level of understanding, since it affects both the evolution of life and star formation. 
However, gravity does not change chemistry, which is already a big simplification. So it is 
not inconceivable that the value of Newton's constant may in the future receive an anthropic 
explanation. 

In this paper we discussed anthropic approaches to solving the two cosmological constant 
problems (CCPs). The first (old) CCP is the discrepancy between the observed small value of 
Pa and the large values suggested by particle physics models. The second (time coincidence) 
CCP is the puzzling coincidence between the epoch of galaxy formation t G and the epoch of 
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A-domination t\. While it is conceivable that the old CCP can be resolved by fundamental 



physics [|42] , |43i |7l,f44|1 , we have argued that the time coincidence problem calls for an anthropic 
explanation. 

We first considered models with a discrete spectrum of p\. These include Abbott's scalar 
field model with a "washboard" potential as well as models in which pa can change 
through brane nucleation accompanied by a change of the four-form field F p5|,|5|j6|.|9[1 . Such 
models can solve both CCPs, provided that (i) the separation between the discrete values 
of pa is e ~ Pmo, where pmo is the present matter density, (ii) the probability distribution 
for pa at the end of inflation is nearly flat, P*(pa) ~ const, and (Hi) the brane nucleation 
rate is sufficiently low, so that the present vacuum energy does not drop significantly in less 
than a Hubble time. We discussed the cosmology of this class of models, the calculation of 
the prior distribution V*(pa), and the observational constraints on the model parameters. 

The required values of the "level separation" e may appear uncomfortably small, but 
Feng, March- Russell, Sethi and Wilczek (FMSW) || have argued that they can naturally 
arise due to non-perturbative effects in M-theory. In M-theory-related models, the brane 

1 10 

tension a is related to e through e ~ ap bare /M p and should also be very small. Our analysis 
shows that in such models the conditions (i)-(iii) cannot be satisfied without fine-tuning of 
the parameters. 

It was conjectured by Weinberg [12] that the condition (ii) of a flat a priori distribution 
for A would automatically be satisfied in any particle physics model where the cosmological 
constant is a random variable. In Ref. |4j] we showed that this conjecture is not always 
satisfied in models where the role of the cosmological constant is played by a slowly varying 
field. Here, we have shown that the conjecture is generically not satisfied in four-form models 
either. In fact, this condition has to be enforced in order to fit observations. This, in turn, 
places severe constraints on the model parameters. Hence, in trying to solve the cosmological 
constant problems by anthropic means, the flat a priori distribution for A cannot be taken 
for granted and the problem of calculating V* has to be addressed. 

Bousso and Polchinsky || have studied models with multiple four-form fields Fi and 
found that the spectrum of the allowed values of pa can be sufficiently dense even for large 
brane tensions. However, in this case the vacua with nearby values of pa have very different 
values of F^ and a flat probability distribution required in (ii) is rather unlikely. Moreover, 
the low-energy physics in such vacua is likely to be different, and it appears that the anthropic 



approach to solving the CCPs cannot be applied in this case |26[ . 

For models unrelated to M-theory, a and e are generally unrelated, and values consistent 
with the constraints (i)-(iii) can easily be found. However, if one gives up the M-theory 
connection, then the FMSW argument cannot be used, and one has to seek an alternative 
explanation for the tiny value of e. Alternatively, one might seek modifications of the FMSW 
model that could relax the relation between e and a. 

All of the earlier discussions of the cosmology of discrete A models encountered the 
"empty universe problem" |31 , 19 , 25| ,|^[6||||26[] . In order to make the present vacuum suffi- 



ciently stable, the brane nucleation has to be strongly suppressed. One then finds that the 
time it takes the universe to evolve from some initial high value of pa to the present low value 
is much greater than the present Hubble time. This suggests that by the time the process 
is complete, any matter that the universe initially had may get diluted to an extremely low 
density, so that one would end up with an empty universe dominated by the cosmological 
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constant. 

We have argued that the empty universe problem disappears when the eternal nature 
of inflation is taken into account. During inflation, brane nucleations leading to higher 
and lower values of p\ have nearly equal probabilities. As a result, the values of are 
randomized, with different parts of the universe thermalizing with different values. The 
resulting probability distribution V*(pa) can be calculated using the stochastic formalism 
we developed in Section IV. The slow rate of brane nucleation is not a problem in eternal 
inflation, since an unlimited amount of time is available. 

FMSW suggested an interesting possibility that nucleation of multiple branes could be 
enhanced by a large degeneracy factor due to the light fields living on the branes. If true, 
this could significantly modify the brane model cosmology. In Section V we studied multiple 
brane nucleation in some detail and found that the pre-exponential factor in the brane 
nucleation rate can both enhance and suppress multiple brane nucleation, depending on 
the field content of the branes. We also concluded that models in which multiple brane 
nucleation dominates can be ruled out, because in such models there is nothing to prevent 
our present vacuum from tunneling down to deep anti-de Sitter space. 

We also discussed models with a continuous spectrum of p\, in which the role of the 
cosmological constant is played by the potential V(<f>) of a scalar field <j)(x). The potential 
has to be very flat, so that its value does not significantly evolve on the present Hubble time 
scale. The values of the field <ft are randomized by quantum fluctuations during inflation, 
and models can easily be constructed in which the resulting probability distribution for 
V(4>) is nearly flat in the range of interest, thus solving both CCPs. The challenge here is to 
justify the very flat potentials required in this class of models. Possibilities include a pseudo- 
Goldstone field which acquires a potential through instanton effects [Q, a large running of the 
field renormalization WE\ , and a non- minimal kinetic term with an exponential ^-dependence 
|| . We have pointed out some difficulties of the instanton approach. 

We thus see that both discrete and continuous A models could in principle solve both of 
the CCPs. However, none of the models that have been suggested so far appears particularly 
well motivated or natural. 

An alternative approach is to assume that the old CCP can be solved within the fun- 
damental theory. The cosmological constant is then truly a constant and is given by an 
expression like pa ~ M^Mp 2 , as in [[?]]. At the same time, the amplitude of density fluctua- 
tions Q could be a random variable, so that the epoch of galaxy formation to is different in 
different parts of the universe. We have shown in Section VIII that, for a wide class of prior 
distributions most of the galaxies will be in regions where t G ~ £a, thus explaining 

the cosmic time coincidence. It would be interesting to extend this analysis and calculate 
the distribution V*(Q) for some models with a variable Q. One would then have some idea 
of how naturally the distributions of the required type can be obtained. 
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NOTES ADDED 



1- After this paper was submitted for publication, a revised version of Ref. || has ap- 
peared. There, it is pointed out that the relation (21) between the tension a and the charge 
density q of the brane does not hold for branes wrapping on degenerating cycles. Instead, the 
tension is suppressed by an exponential factor relative to the charge. We note two potential 



problems with this picture. First, as it was argued in Ref. fl45 |, the brane charge and tension 



appear to be unprotected against renormalization below the supersymmetry breaking scale. 
Such renormalization would make the brane charge q unacceptably large. Second, if for some 
reason the brane parameters do not get renormalized, then, in order to satisfy the anthropic 
constraint (22) on q, the brane tension has to be exceedingly small. The instanton action 
(9) would then be small and brane nucleation would be completely unsuppressed. 

2- A new approach to explaining very flat scalar potentials and branes with a very small 
four- form charge has been suggested in Ref. |4"o| , where these features are attributed to a 
spontaneously broken discrete symmetry. 
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